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Abstract 



The excitation curve for the r"*" r~ production in electron positron annihilation near 
the threshold is revisited with the aim of updating and extending a previous work. We 
find that the corrections of the relative magnitude 0{a) near the threshold are signifi- 
cantly contributed by the radiative modification of the Coulomb interaction between the 
r leptons. The interpolation between the Coulomb effects at the threshold and the rela- 
tivistic effects at higher energies is considered and the resulting formula is argued to have 
relative accuracy up to, but not including, terms of the order of at any velocity of the 
T leptons. 



The dramatic improvement in the accuracy of the measured mass of the r lepton, 
achieved in the BES experimentfl, has demonstrated one of the advantages of studying 
the production of the r"*" r~ pair in e+ e~ annihilation in the immediate vicinity of the 
threshold. In view of the special kinematical and background advantages of the threshold 
regionS one may believe that experiments at the r+ threshold will be continued at the 
existing machine and, hopefully, at the much anticipated tau-charm factoryS. In view 
of the probable improvement in the precision of the experimental data, it is worthwhile 
to have sufficiently accurate theoretical description of the e~ r+ r~ cross section 
at the threshold and beyond. In addition, this process is of a theoretical interest of its 
own, since it provides a 'test ground' for QED almost as good as the muon or electron 
magnetic anomaly. 

The well-known lowest-order QED formula for the cross section 

a,{e^e- ^r^r-) = ^vi3-v') , (1) 
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where v = yl — 4m^ /s is the velocity of each of the r leptons in the cm. frame, is subject 
to corrections, arising from various sources: 

i - radiation from the initial electron and positron, 

ii - vacuum polarization in the time-like photon, 

m- corrections to the spectral density p(g^) = — | l]x(0|j^(— 5')|-^)(-^|j/i(q')|0) of the 
electromagnetic current = {f'-f^r) of the tau leptons. The latter corrections generally 
include both the electromagnetic interaction between the r leptons and the final state 
radiation. 

In higher orders there appears interference between the effects i - in, however its 
relative magnitude is suppressed by at least the factor a^, which makes it beyond the 
scope of the present paper. 

Because of the radiative effects the actually measured cross section is given by the 
formula 

/■W 

a{W)= r{W,w)\l~U{w)\-'^a{w)dw , (2) 

where W = ^/s is the total energy in the cm. system. The weight function r{W,w) 
describes the radiation from the initial state@ and |1 — n(w)|~^ is the vacuum polarization 
factori. These two effects are standard and are routinely accounted for in the data 
analyses, while the dynamics of the final state is encoded in the cross section a{w) = 



1 



In this letter we concentrate on the behavior of (y{w) in the threshold region, where 
the velocity v of the produced r's is a small parameter. In this region the QED effects 
in the spectral density p(if^) are determined by the Coulomb interaction between the 
r leptons at distances r ~ (mT-f)~^, the parameter for which effects is given by vra/f, 
while the relativistic effects and the radiation of transversal photons by the r leptons 
are suppressed by at least the factor v^. In the non- relativistic region one can effectively 
separate the dynamics at short distances r ~ Tn~^ , at which the leptons are produced by 
the electromagnetic current and that at the distances r ~ (mT-f)~^, which determine the 
Coulomb effects, by replacing the current by a local non-relativistic operator. In the 
lowest order in the loop corrections, arising at short distances, this replacement can be 
written as (c.f. Ref. 0) 

a{w) = — ImG(0,0;m^t;2) , (3) 

where (^(x, y; E) is the non-relativistic Green's function at energy E = w — 2m^. For free 
particles the Green's function is given by 

m An |x — y| 

so that ImGolO, 0; m,- f ^) = im?v/{An) and the non-relativistic limit of the 'bare' cross 
section (|l|) is reproduced. The solution for the Green's function in the Coulomb poten- 
tial V{r) = —a/r is also well known and at the origin it reads as ImGc(0, 0; f ^) = 
Fc ImGo(0, 0; t>^) with the Coulomb factor Fc given by 

(5) 



1 — exp(— vra/f ) 

The behavior of this factor at v ^ makes the cross section at v = finiteS and of 
a measurable magnitude: 

a(v = 0) = (1 + 0(a)) = 0.236 nb (1 + 0(a)) . (6) 

The fact that the spectral density of the current at small velocity is determined by the 
Coulomb interaction calls for the idea that the radiative corrections to this interaction 
should show up at the level of 0{a) rather than the usual O(a^). Since the higher order 
effects, involving exchange and radiation of transversal photons at large distances, are 
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suppressed by the factor t>^, the only radiative effect at these distances is due to running 
of the couphng a, which is also known as the Uehling-Serber radiative correction to the 
potential (see e.g. in the textbook j^). For a preliminary estimate of the magnitude of 
this effect one can replace a in the Coulomb factor (j^) by the effective coupling at the 
momentum p = m^f : 

a ^ a 1 + — m . (7) 



One can readily see that the ratio rriT-v/me is large even for small v, so that the effect can 
be of the order of one percent. This estimate however requires some elaboration, which 
is the main subject of the present work. Namely, non-logarithmic terms can be essential 
and also the estimate of the characteristic distances, essential for the Coulomb factor, is 
modified by the Coulomb interaction at f ~ 0{a). To find the effect quantitatively we 
use the full form of the Uehling-Serber correction to the Coulomb potential due to the 
electron-positron vacuum polarization^: 



and consider this correction as a perturbation, thus finding the modification of the wave 
function at the origin in the form 



6iIj{0) = - j Gc{0,r;my)6V{r)^c{r)d^r , (9) 
where ipcif) is the S-wave wave function at energy E = itLt-v'^ in the Coulomb field —a/r: 

ilj^{r) = C e-'P' iFiil + i\2,2ipr) (10) 

with p = rrirV, A = mr<y/{2p) = a/{2v) and C being normalization constant. Using the 
representation of the Coulomb Green's function in the form 

GMr;pVm^) = -^l^e^^'- J\'^^^' (^)^' dt , (11) 

we find that the corrected expression for the square of the wave function at the origin can 
be written as: 

\Mo) + sm\' = \Mo)\'a + h) , (12) 

where the correction h after integration over r in eq.(0) is given by 
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2a 
3^ 



-2A Im / dt \ dx{ —— ] ^ ' ' ' ' ^ 

Ji 



(13) 

with z = me/rrir- 

The double integral in the latter equation was calculated numerically^. The plot of 
the magnitude of the relative correction /i as a function of velocity v is shown in Fig.l 
along with the result of the simple logarithmic estimate from the equations (^ and (0). 
One can see that this estimate becomes invalid in the region v < 0{a), while it fits the 
exact result quite well at higher velocity. Also in Fig.l is shown the result of calculation 
of the contribution to h of the muon vacuum polarization, which, as expected, is quite 
small. We also find that in the limit f — > and me/mr the quantity in the square 
brackets in eq.(|l^) is equal to 4, so that h{v = 0,me = 0) = 8a/(37r). 

In order to have a full description of the 0{a) correction to the excitation of the r"*" r~ 
channel near the threshold one should also take into account the finite renormalization 
of the electromagnetic current, which comes from short distances r ~ m~^, which is well 
known 01 to result in the overall factor (1 — Aa/rr), which, unlike the term h, does not 
depend on the velocity in the threshold region, i.e. unless the terms 0(f ^) are taken into 
account. 

Summarizing this discussion we find that the cross section for production of non- 
relativistic r leptons with velocity v including the terms of the relative magnitude 0{a) 
is given by 

^ = ^1 f + ^ ' 

2m:; 1 — exp(—7r a/t>j V tt / 

where h is determined by the running of the QED coupling a and the dominant contribu- 
tion to it due to the electron-positron vacuum polarization is given by eq.(pr3D and which 
is plotted on Fig.l. 

At the present accuracy of experimental measurement of rrir, which is about 0.3 
MeVfl, the 0{a) corrections discussed here can be completely ignored. A simple analysis 
shows that these terms become essential for measurements with accuracy of the order 
of 10^^ MeV, at which level one should also take into account the contribution of the 



^In particular we found that the numerical integration routines built into Mathematica provide a 
reliable, though not the fastest, platform for this calculation. 
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Coulomb bound stateso. For completeness we cite here that the binding energy of the 
n-th level is 



\E\ = , ^ ^ 5 — (15) 



and its e~^e decay rate is given by 

rriT-a^ 6.1 eV 



« 3 - 3 • (16) 

Equation (0) is applicable in the non-relativistic limit, i.e. the relative magnitude of 
the corrections to it is determined by v"^. On the other hand the exact in v formula of 
the first order in a is known (see e.g. in Schwinger's textbook [^): a{e^e^ t^t^) = 
ao (l + ^S{v)^, where ctq is the bare cross section, eq.(|l]), and 

Siv) = + [t + 1^ i"^) 1^ (t^^) + 2 {^) + 2 Li, (1±H) - 

2 Li2 {^) - 4 Li^iv) + U,{v')] + [^(1 + v') -3v+ ^j^] In (i±f) + 

(17) 

with Li2(x) = — /g^ln(l — t) dt/t = I]$^Li ^^"'/n^. The full relativistic formula of the first 
order in a develops at low velocity the inaccuracy of the order of (a/f )^, which matches 
the inaccuracy O(f^) of the non-relativistic formula (|14|) at v"^ ~ a, where the relative 
magnitude of the error in either of these formulas is O(f^) ~ 0{a). One can however write 
an expression 01 , which interpolates between the non-relativistic formula (p^, summing 
all the terms of the form (a/v)'' and a (a/f)^, and the full formula of the first order in 
a. The resulting interpolating formula can be argued to have relative accuracy 0{a'^) at 
any velocity]^ 

The interpolating formula has the form 

a{e+e- ^ r+r~) = (1 + h) (l + ^S{v) - ^) (18) 

where the lowest-order cross section (Tq is given by eq.([l|), the Coulomb factor is in 

the equation (|^), the correction term h, given by eq.(|l3]), accounts for running of a in 

the Coulomb terms and the function S{v) is written in eg. ([171). One can notice that in 

eq.(|18D the first Coulomb term 7Ta/{2v) is subtracted from the full relativistic expression 

^We take this opportunity to correct some misprints and other minor flaws of the discussion of the 
interpolating formula in Ref. . These flaws however do not affect the results of the data analysis in the 
BES experiment 01 , which used the formulas from Ref.fl. 
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of the first order in a, since this term is contained in the factor F^. The factor (1 — 4a/7r), 
coming from finite renormahzation of the current at the threshold, which is present in the 
equation (p!^, in the interpolating formula is contained in the factor with S{v). Equation 
(^) sums all terms of the form {a/vY and a {a/v)^ as well as all terms linear in a with 
arbitrary dependence on v. Therefore the corrections to the interpolating formula start 
with terms of the form (a/f)^ with A; > 0, which are at least as small as O(a^) in 
comparison with those included in eq. ([l8|) , independently of the velocity v. Naturally, to 
also take into account the terms of the relative order one needs a full calculation of the 
two-loop correction to the spectral density of the electromagnetic current at arbitrary v 
as well as a calculation of the interference between the mechanisms i - Hi mentioned in 
the beginning of this paper. 

This work was supported, in part, by the DOE grant DOE-AC02-83ER40105. 
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Figure caption 

Figure 1. The relative correction h (eq.(|TEp) to the cross section of r+r^ production 
near the threshold, arising from modification of the Coulomb potential by the polarization 
of the electron-positron vacuum (solid) and of the muon one (dashed) vs. the velocity of 
the r leptons in the cm. system. Also is shown (dotted) the logarithmic approximation 
for the electron-positron contribution to h according to eq.(^. 
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